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Abstract 

In this survey, different aspects of the theory of orthogonal polynomials 
of one (real or complex) variable are reviewed. Orthogonal polynomials 
on the unit circle are not discussed. 
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1 Introduction 

The theory of orthogonal polynomials can be divided into two main but only 
loosely related parts. The two parts have many things in common, and the 
division line is quite blurred, it is more or less along algebra vs. analysis. One 
of the parts is the algebraic aspect of the theory, which has close connections 
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with special functions, combinatorics and algebra, and it is mainly devoted 
to concrete orthogonal systems or hierarchies of systems such as the Jacobi, 
Hahn, Askey-Wilson, ... polynomials. All the discrete polynomials and the q- 
analogues of classical ones belong to this theory. We will not treat this part; 
the interested reader can consult the three recent excellent monographs [39] by 
M. E. H. Ismail, [28] by W. Gautschi and [6] by G. E. Andrews, R. Askey and 
R. Roy. Much of the present state of the theory of orthogonal polynomials of 
several variables lies also close to this algebraic part of the theory. To discuss 
them would take us too far from our main direction; rather we refer the reader 
to the recent book [24] by C. F. Dunkl and Y. Xu. 

The other part is the analytical aspect of the theory. Its methods are ana¬ 
lytical, and it deals with questions that are typical in analysis, or questions that 
have emerged in and related to other parts of mathematical analysis. General 
properties fill a smaller part of the analytic theory, and the greater part falls 
into two main and extremely rich branches: orthogonal polynomials on the real 
line and on the circle. The richness is due to some special features of the real 
line and the circle. Glassical real orthogonal polynomials, sometimes in other 
forms like continued fractions, can be traced back to the I8th century, but their 
rapid development occurred in the 19th and early 20th century. Orthogonal 
polynomials on the unit circle are much younger, and their existence is largely 
due to Szegd and Geronimus in the first half of the 20th century. B. Simon’s 
recent treatise [80, 81] summarizes and greatly extends what has happened since 
then. 

The connection of orthogonal polynomials with other branches of mathe¬ 
matics is truly impressive. Without even trying to be complete, we mention 
continued fractions, operator theory (Jacobi operators), moment problems, an¬ 
alytic functions (Bieberbach’s conjecture), interpolation. Fade approximation, 
quadrature, approximation theory, numerical analysis, electrostatics, statistical 
quantum mechanics, special functions, number theory (irrationality and tran¬ 
scendence), graph theory (matching numbers), combinatorics, random matrices, 
stochastic processes (birth and death processes; prediction theory), data sorting 
and compression. Radon transform and computer tomography. 

This work is a survey on orthogonal polynomials that do not lie on the unit 
circle. Orthogonal polynomials on the unit circle—both the classical theory and 
recent contributions—will be hopefully dealt with in a companion article. 

This work is meant for non-experts, and it therefore contains introductory 
materials. We have tried to list most of the actively researched fields not directly 
connected with orthogonal polynomials on the unit circle, but because of space 
limitation we have only one or two pages on areas where dozens of papers 
and several books had been published. As a result, our account is necessarily 
incomplete. Also, the author’s personal taste and interest is reflected in the 
survey, and the omission of a particular direction or a set of results reflects in 
no way on the importance or quality of the omitted works. 

For further backgound on orthogonal polynomials, the reader can consult 
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the books Szego [91], Simon [80]-[81], Freud [27], Geronimus [34], Gautschi [28], 
Chicara [18], Ismail [39]. 

This is a largely extended version of the first part of the article Golinskii- 
Totik [36]. 

Acknowledgement. Research was supported by NSF grant DMS-040650 and 
OTKA T049448, TS44782, and was carried out within the Analysis Research 
Group of the Hungarian Academy of Sciences. The author participates in the 
project INTAS 03-51-6637 that supported a visit by Leonid Golinskii to Szeged, 
during which the outline of this paper was laid out. 


2 Orthogonal polynomials 

Orthogonal polynomials with respect to measures 

Let /i be a positive Borel measure on the complex plane, with an infinite number 
of points in its support, for which 

J \z\'^dfj,{z) < oo 

for all m > 0. There are unique polynomials 

Pniz) = Pnip-, Z) = H-, > 0, 71 = 0, 1, . . . 

that form an orthonormal system in i.e. 

f _, f 0 if m 7 ^ n 

J I 11 m = n. 

These p„’s are called the orthonormal polynomials corresponding to p. Kn is 
the leading coefficient, and Pn{z)/Kn = 2 ;" -I- • • • is called the monic orthog¬ 
onal polynomial. The leading coefficients play a special and important role in 
the theory, many properties depend on their behavior. When dp(x) = w(x)dx 
on some interval, say, then we talk about orthogonal polynomials with respect 
to the weight function w. 

The pnS can be easily generated: all we have to do is to make sure that 

f Pn{ ) fc = 0,1,..., 77 — 1, (2.1) 

J 

which is an 77 X 77 system of equations for the non-leading coefficients of Pn{z) j Kn 
with matrix where 




z'^z^dp(z) 
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are the complex moments of /x. This matrix is nonsingular: if some linear 
combination with coefficients Cq, . •., c„_i of the rows is zero, then the polyno¬ 
mial Pn-i{z) = Co -b • • • -b is orthogonal to every , j < n, and hence 

it is orthogonal to itself, i.e., 

J \Pn-l\'^d^i = J Pn-lPn-ldfi = 0, 


which implies Pn-i(z) = 0. Thus, co = ••• = c„_i = 0, which shows the 
nonsingularity of {crij). Therefore, the system (2.1) has a unique solution for 
the non-leading coefficients of Pn{z)/Kn (note that the leading coefficient is 1), 
and finally comes from normalization. 

In particular, the complex moments already determine the polynomials. In 
terms of them one can write up explicit determinant formulae: 


Pn{z) = 




Co.o 

CTO.l 

Cr0,n-1 

I 

Cl.o 

CTl,l 


z 

^n-1,0 

CTn-l,! 

' ' ^n—l,n—l 

ZU-I 

^n,0 

C'n,l 

^n,n — l 

z” 


where 


Pn. — \o'i. 


'^,3 \i,j=0 

are the so called Gram determinants. 

Note that if /x is supported on the real line then 


( 2 . 2 ) 


(2.3) 




J x*+^(i/x(x) =: ai+j, 


so Dn = \ai+j\f j^Q is a Hankel determinant, while if /x is supported on the unit 
circle then 


^i,3 


z" ^dp{z) =: 


so Dn = \f3i-j\2j^Q is a Toeplitz determinant. In these two important cases the 
orthogonal polynomials have many special properties that are missing in the 
general theory. 


The Riemann Hilbert approach 

Let /X be supported on the real line, and suppose that it is of the form c?/x(t) = 
w(t)dt with some smooth function w. A new approach to generating orthogonal 
polynomials that has turned out to be of great importance was given in the early 
1990’s by Fokas, Its and Kitaev [26]. Consider 2x2 matrices 

Y{z) = 


Yn{z) Yi 2 {z) 
Y2iiz) Y22 {z) 
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where the Yij are analytic functions on C \ R, and solve for such matrices the 
following matrix-valued Riemann—Hilbert problem: 

1. for all a; € R 

r-(x) = y-(x)(J ”«) 

where Y~^, resp. Y~, is the limit of Y(z) as z tends to x from the upper, resp. 
lower half plane, and 

yw = (/ + o(l))(; 

at infinity, where I denotes the identity matrix. 

Thus, one is looking for 4 functions Yu,..., Y 22 analytic on C\R, such that if 
denote the boundary limits of these functions at x G R from the upper, 
resp. lower half plane, then 

^11 (^)) ^21 (2^) = ^21(2^) ( 2 - 4 ) 

and 


^ 12 ( 2 ;) = 4^n(a;)w(x) -b Y{^{x), Y^^{x) = Y^^{x)w{x) + Y2^(x). (2.5) 


These connect the functions on the upper and lower half planes only very mildly, 
and what puts the problem into a rigid form is the second condition, namely it 
is required that for large z uniformly on the plane we have 

Yn(z) = z" + 0(|zr-i), Y2i(z) = 0(|zri) (2.6) 

and 

Yi 2 (z) = 0(z-”-i), Y 22 (z)=z-" + 0(|z|-"-1). (2.7) 

It can be shown that there is a unique solution Y(z). The relevance of this to 
orthogonal polynomials is that the entry Yii(z) is precisely the monic polyno¬ 
mial Pn(/r, z)lKn- Indeed, (2.4) shows that Yu and Y 12 are analytic everywhere, 
and if an entire function is Odzl™) as z ^ 00 , then it is a polynomial of degree 
at most m. Thus, we get from (2.6) that Yii(z) = z"-|- - • • is a monic polynomial 
of degree n, and Y 2 i(z) is a polynomial of degree at most n — 1. The relation 
(2.7) gives that the integral of z*Yi 2 (z) over the circle |z| = R is 0(i?^“”) for 
dX\ k < n and hence it tends to 0 as i? —> 00 . By analyticity, the integral over 
the upper part of the circle can be deformed into an integral from R to —R on 
the upper part of R, i.e., into 




x^Y^{x)dx, 


and similarly the integral over the lower part of the circle can be deformed into 
an integral from —R to R on the lower part of R, i.e., into 

l-R 

/ x^Y^2{x)dx. 

J-R 
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The first relation in (2.5) implies 

x^Y{^{x) — x^Y^{x) = x^Yii{x)w{x), 


therefore for fc = 0, 1 we have 

rR 


[ x'^Yn{x)w{x)dx = 0{R^-^) = 0{R-^) 0, 

J-R 


which implies 



x’^Yii{x)w{x)dx = 0 . 


Thus, 111 is indeed the monic n-th orthogonal polynomial with respect to w. 

The other entries can also be explicitly written in terms of the orthogonal 
polynomials and Pn-i- Y 21 is a constant multiple of Pn-i, 


Yi2{z) 


1 

2iTTKn 


Pn{x)w{x) 

- dx 

x — z 


is the Cauchy transform of Pn{x)w{x) / Kn, and Y 22 is the Cauchy transform of 
y 2 i (= const-pn-i). Furthermore, and the recurrence coefficients a„,5„ (see 
Section 4) can be expressed in terms of the entries of Yi , where Yi is the matrix 
defined by 

For details on this Riemann-Hilbert approach, see Deift [20]. 


Orthogonal polynomials with respect to inner products 

Sometimes one talks about orthogonal polynomials with respect to an inner 
product (•,•) which is defined on some linear space containing all polynomials, 
and orthogonality means {pmPm) = 0 for m ^ n. In this case the aforemen¬ 
tioned orthogonalization process can be used, and with the 

determinantal formula (2.2) is still valid. 

Sometimes one has an (•,•) with the standard inner product properties, 
except that positive definiteness may not hold (as an example consider non- 
Hermitian orthogonality from Section 14). Then the orthogonalization process 
and the determinantal formulae can still be used provided the Gram determi¬ 
nants (2.3) are different from zero. If this is not so, then we write 

Pn{z) = 7n2” -h H-, 

and make sure that Pn is orthogonal to all powers , 0 < k < n, i.e., solve the 
homogeneous system of equations 

n 

= 0 , 

1=0 


fc = 0,..., n — 1, 
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for 7 o, 7 i,..., 7 „. Since the number of unknowns is bigger than the number of 
equations, there is always a non-trivial solution, which gives rise to non-trivial 
orthogonal polynomials. However, then we cannot assert any more 7 „ 7 ^ 0, so 
the degree of may be smaller than n, and there may be several choices for 
Pn- Still, in applications where non-Hermitian orthogonality is used, these Pn 
play the role of orthogonal polynomials. 

Varying weights 

In the last 25 years orthogonal polynomials with respect to varying measures 
have played a significant role in several problems, see e.g. the sections on expo¬ 
nential and Freud weights or on random matrices in Section 4. In forming them 
one has a sequence of measures /x„ (generally with some particular behavior), 
and for each n one forms the orthogonal system {pkipn, most cases 

one needs the behavior of Pn^Pn^z) or that of _Pn±fc(/r„,z) with some fixed k. 
We mention three examples. 

The first example is that of Freud weights: W{x) = , A > 0. If one 

substitutes x = n^/^y, then with Pn{y) = Pn{W, x) orthogonality takes the form 

y = 0, n ^ m, 

and it turns out that this is just the right scaling, e.g. the zeros of P„ have 
an asymptotic zero distribution (while those of p„(IF,z) are spreading out to 
infinity). Thus, studying orthogonal polynomials with respect to Freud weights 
W is equivalent to studying orthogonal polynomials with respect to the varying 
weights Wn(x) = IF(x)", and actually, working with Wn turns out to be very 
natural. 

For the second and third examples see multipoint Fade approximation and 
random matrix theory in Section 6 . 

Matrix orthogonal polynomials 

Orthogonality of matrix polynomials (i.e., when the entries of the fixed size ma¬ 
trix are polynomials of degree n = 0 , 1 ,... and orthogonality is with respect to 
a matrix measure) is a very active area which shows extreme richness compared 
to the scalar case. See Section 16 for a short discussion. 


3 The extremal problem 

One of the most useful tools in the study of orthogonal polynomials is the fact 
that they solve the following extremal problem: minimize the norm for 

all monic polynomials Pn{z) = + ■ ■ ■ of degree n. The minimum turns out 
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to be 1 /k^, i.e., the n-th monic orthogonal polynomial is the (unique) extremal 
polynomial in 

inf [\Pn\^dfi=^. (3.1) 

P„(z)=z" + -.7 

Indeed, any Pn is a linear combination X]fe=i CkPk with the orthonormal polyno¬ 
mials pfc, and here c„ = I/kh because Pn is monic, i.e., it has leading coefficient 
1. Orthogonality gives 

/ n 

\Pn\^dp = ^ |CfcP, 
fc=0 

from which we can see that this is always > jcnp = 1 /k^, and equality occurs 
if and only if all the other Cfc’s are 0. 

A related extremum problem leads to the so called Christoffel functions 
associated with p. They are defined as 


Xn{p,z)= inf [ \Pn\‘^dp. 

P„{z) = l, deg{P„)<uJ 

If we write again Pn = O? then Pn{z) = 1 means 

n 

'^CkPk{p,z) = 1 , 
k^O 

and hence by Cauchy’s inequality 

n n 

1 < icfen(^ \pk{p,z)\^). 

k^O k^O 

Therefore, 

/ n n 

\Pn\^dfi = Y icfcp > Y1 

1 _n 7,_/-» 






with equality if and only if 


^k — 


Pkip,z) 


ELo \PkiP,zW 
Thus, we have arrived at the formula 


(3.2) 


Xnip,z) ^ = y] |pfc(Ai,z)P (3.3) 

fc=0 

for all z G C for the Christoffel function A„(^, z). 

For example, for measures p lying on the real line it is easy to see from 
this formula that p has a point mass at xq, i.e., /x({xo}) > 0 if and only if 
J2kPk{pTXo)‘^ < oo, and then 

OO 

P{{xo}) = {YPk{P,Xo)^)~^. 
fc=0 
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4 Orthogonal polynomials on the real line 


Let /X be supported on the real line. In this case orthogonalization leads to real 
polynomials (i.e., all the coefficients are real). The most remarkable property 
of this real case is that the p„’s obey a three-term recurrence formula 


XPn{x) = anPn+l{x) + b^Pnix) + On-lPn-lix), 


where 


^r?, — 


^n+1 


> 0, ^ / xpl{x)dp{x) 


(4.1) 


(4.2) 


are the so called recurrence coefficients. Indeed, if we write xpn{x) as a 


linear combination ^kPk(z) with 

Ck ■■= j xpn{x)pk{x)dp{x), 

then all the c^’s for fc < n — 1 vanish by orthogonality: 

Ck = J xpn{x)pk{x)dfj,{x) = J {xpk{x))pn{x)dfj,{x) = 0 


because xpk(x) is a polynomial of degree smaller than n. Comparison of the 
leading coefficients on both sides gives that c„+i = KniKn+i, but since c„+i is 
also the integral of xpn{x)pn+i(x) against /i, we get that 

C„_1 = / xpn{x)pn-l{x)dp{x) = _ 

J Kn 

Finally, c„ is the integral given in (4.2). 

We emphasize that the three-term recurrence is a very special property of 
real orthogonal polynomials, and it is due to the fact that in this case the 
polynomials are real, hence 


J Xpn{x)pk{x)dp{x) 


J pn{x){xpk{x))dp{x) = 0 


for fc < n — 1. In the non-real case the two sides here are totally different. 
The three-term recurrence is missing in the general case, and it is replaced by 
a different recurrence for polynomials on the circle. 

Conversely, any system of polynomials satisfying (4.1) with real a„ > 0, 
is an orthonormal system with respect to a (not necessarily unique) measure on 
the real line (Favard’s theorem). The unicity of the measure in question is the 
same as the determinacy of the moment problem, which in turn is again closely 
related to the behavior of orthogonal polynomials; see Section 6. 

In the real case the zeros of are real and simple and the zeros of 

and Pn+i interlace, i.e., in between any two zeros of Pn+i there is a zero of 
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Pn- In fact, Pn must have n sign changes, for if it had only m < n, say at the 
points j/i, G R, then it could not be orthogonal to the polynomial q{x) = 

Y[jLii^~yj) of degree m < n, for then q{x)pn{x) would be of constant sign. Let 
now Xn < Xn-i < ■ ■ ■ < xi he the zeros ofp„, and suppose that we already know 
that the zeros of Pn andp„_i interlace, which implies sign(p„_i(xfc)) = (—1)^“^. 
If we substitute Xk into the recurrence formula (4.1) then a„ > 0 gives that 
Pn+i(xk) and Pn-i{xk) are of opposite signs at Xk, i.e., sign(p„ +i{xk)) = i-f)’", 
and this gives that the zeros of and Pn+i also interlace. Thus, the interlacing 
property follows by induction. 

The three-term recurrence implies for the so called reproducing kernel 
the Christoffel-Darboux formula 


'^Pk{x)pk{t) 


Kn Pn+l{x)Pn{t) - Pn{x)pn+l{t) 
X t 


(4.3) 


Indeed, use the recurrence formula for Pn+i on the right and a„ = 
then induction gives (4.3). The special case 


Xn{p,x) ^ =^Pk{xf = {p'n+l{x)Pn{x) - p'n{x)pn+l{x)) (4.4) 

is worth mentioning. 

The starting values of the recurrence (4.1) are p_i = 0, po = (m(C)))~^^^. If 
one starts from g_i = —1, go = 0 and uses the same recurrence (with a_i = 1) 


xqn{x) = anqn+i{x) + bnqn{x) + an-iqn-l{x), (4.5) 


then g„ is of degree n — 1, and by Favard’s theorem the different g„’s are or¬ 
thogonal with respect to some measure. The g„’s are called orthogonal poly¬ 
nomials of the second kind (sometimes for p„ we say that they are of the 
first kind). They can also be written in the form 

J Z - X 


5 Classical orthogonal polynomials 

These are 

• Jacobi polynomials P^’^\ a,/3 > —1, orthogonal with respect to the 
weight (1 — x)°‘{l + x)^ on [—1,1], 

• Laguerre polynomials a > —1, with orthogonality weight x°‘e~^ 
on [0, oo), 
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Hermite polynomials orthogonal with respect to e on the real 
line (— 00 , 00 ). 


In the literature various normalizations are used for them. 

They are very special, for they possess many properties that no other or¬ 
thogonal polynomial system does. In particular, 

• they have derivatives which form again an orthogonal polynomial system, 
e.g. the derivative of is a constant multiple of 


= i(n + a + /3 + 


they all possess a Rodrigue’s type formula 

dnw[x) ax” 

where w is the weight function and ct is a polynomial that is independent 
of n, for example, 


1 (T 


LM(x)=e"x-“-—(e-"x”+“), 

• they satisfy a differential-difference relation of the form 

7r{x)P^{x) = {anX + f3n)Pn{x) + JnPn-l{x), 


e.g. 

xiLi^^Yix) = nL\^\x) -{n + a)L|f_\(x), 

• they satisfy a non-linear equation of the form 

(7(x) {^Pji{x^Pri — \{x^^ — (^CXfiX Pn^Pni^X^Pn—l^x') 

+lnPY{x) + SnP^-lix), 

with some constants an, Pn,'Jn,Sn, and a a polynomial of degree at most 
2, e.g. 

{Hn{x)Hn-i{x))' = 2xHn{x)Hn-i{x) - H^x) + 2nH^_j^{x). 


Every one of these properties has a converse, namely if a system of orthogonal 
polynomials possesses any of these properties, then it is (up to a change of 
variables) one of the classical systems, see Al-Salam [3]. See also Bochner’s 
result in the next section claiming that the classical orthogonal polynomials are 
essentially the only polynomial (not just orthogonal polynomial) systems that 
satisfy a certain second order differential equation. 

Classical orthogonal polynomials are also special in the sense that they pos¬ 
sess a relatively simple 
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• second order differential equation, e.g. 

xy" + (a + 1 — x)y' + ny = 0 


for 

• generating function, e.g. 


E Hn{x) . 2 \ 

-;— w = exp(2a:w — w ), 

n! 


• integral representation, e.g. 

(1 _ ^)«(i+y (1 - t)"+“(i+- x)—^dt 

over an appropriate contour, 

and these are powerful tools to study their behavior. 

For all these results see Szegd [91]. 


6 Where do orthogonal polynomials come from? 

In this section we mention a few selected areas where orthogonal polynomials 
naturally arise. 


Continued fractions 


Continued fractions played an extremely important role in the development of 
several branches of mathematics, but their significance has unjustly diminished 
in modern mathematics. A continued fraction is of the form 

Bi 


and its n-th convergent is 

^ 

Rn Al 


n= 1,2,.... 


^2 + ' 


Bn 


The value of the continued fraction is the limit of its convergents. The denom¬ 
inators and numerators of the convergents satisfy the three-term recurrence 
relations 


Sn — T BnSn—2i 


Rq = 1 , R—i = 0 
5o = 0, = 1, 
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which immediately connects continued fractions with three-term recurrences 
and hence with orthogonal polynomials. 

But the connection is deeper than just this formal observation. Many el¬ 
ementary functions (like z — \/z'^ — 1) have a continued fraction development 
where the B^s are constants while the An's are linear functions, in which case 
the convergents are ratios of some orthogonal polynomials of the second and 
hrst kind. An important example is that of Cauchy transforms of measures fi 
with compact support on the real line (so called Markov functions): 

/(,)= /AW (6.1) 

J X - z z z^ 

The coefficients aj in the development of (6.1) are the moments 

aj = ( x^diJ,{x), j = 0,1,... 


of the measure y,. The continued fraction development 


fiz) 


Bi 


z — Ai -\- 


B2 

Z — A 2 + --- 


of / at infinity converges locally uniformly outside the smallest interval that 
contains the support of /i (A. Markov’s theorem). 

As has been mentioned, the numerators S'„(z) and the denominators i?„(z) 
of the n-th convergents 


Sn{z) 

Rniz) 


Bi 


z — Ai 


_ B2_ 

Z-A 2 + --- 



n = 1,2,... 


satisfy the recurrence relations 


Rn{z) = {Z - An)Rn-l{z) + BnRn- 2 {z), Ro=l,R-l=0 (6.2) 
Sn{z) = {Z - An)Sn-l{z) + BnSn- 2 {z), So = 0, S-i = l. 


These are precisely the recurrence formulae for the monic orthogonal polyno¬ 
mials of the hrst and second kind with respect to /i, hence the n-th convergent 
is cg„(z)/p„(z) with c = ^(C)^/^. 

See Szegd [91, pp. 54-57] as well as Kruschev [42] and the numerous refer¬ 
ences there. 


Fade approximation and rational interpolation 

One easily gets from the recurrence relations (6.2) that 


Rm{z) Sni+l{z') 
Rmiz) Rm+liz'j 


(- 1 )" 


BiB 2 - ■ • Bn+l 
Rn{z)Rn+l{z) 
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and summation of these for m = n,n + 1,... yields that 

2n 


Sniz) -ak , 




2 : 


i.e., with the preceding notation the rational function 
Sn{z)/Rniz) = Cqniz)/pn{z) with C = 


of numerator degree n — 1 and of denominator degree n interpolates f{z) at 
infinity to order 2n. This is the analogue (called [n — 1/n] Fade approxi- 
mant) of the n-th Taylor polynomial (which interpolates the function to order 
n) for rational functions. The advantage of Fade approximation over Taylor 
polynomials lies in the fact that the poles of Fade approximants may imitate 
the singularities of the function in question, while Taylor polynomials are good 
only up to the first singularity. The error in [n — 1/n] Fade approximation has 
the form 


Pn[z) pi[z) 



X — z 


dfj,{x). 


Orthogonal polynomials appear in more general rational interpolation (called 
multipoint Fade approximation) to Markov functions, see e.g. Stahl-Totik 
[87, Sec. 6.1]. For every n select a set An = {xo.n,..., X 2 n,n} of 2n+l interpola¬ 
tion points from C \ / where I is the smallest interval that contains the support 
of p. The points need not be distinct, but we assume that An is symmetric with 
respect to the real line. Fut 


2n 

^n{z) := n (^ - Xjn)- 
i=o 

n ^00 


The degree of tOn is equal the number of finite points in A„. By r„(z) = 
Un{z)/Qn{z) we denote a rational function of numerator and denominator de¬ 
gree at most n that interpolates the function / at the 2n -I- 1 points of the set 
An = {xo^n, • ■ •, X 2 n,n} in the sense that 


= o( 2 -( 2 n+i)) as Jzl^oo; 

UJn{z) 

the expression on the left is bounded at every finite point of An , and at infinity 
it has the indicated behavior. Now for Markov functions this rational inter- 
polant uniquely exists, Qn is the n-th orthogonal polynomial with respect to 
the varying weight dp{x)/uin{x), and the remainder term of the interpolation 
has the representation 


Xln{z) 

Ql{z) 


Qljx) 

U!nix){x - Z) 


dp{x) 


if -rnif, An;-)){z) 
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for all z outside the support of y,. Thus, the rate of convergence of the ra¬ 
tional interpolants is intimately connected with the behavior of the orthogonal 
polynomials with respect to the varying weight dyL{x)/uin{x). 

Moment problem 

The moments of a measure /x, /i(C) = 1, supported on the real line, are 

an = J x^diJ,{x), n = 0, 1,.... 

The Hamburger moment problem is to determine if a sequence {a„} (with nor¬ 
malization Qfo = 1) of real numbers is the moment sequence of a measure with 
inhnite support, and if this measure is unique (the Stieltjes moment problem 
asks the same, but for measures on [0,oo)). The existence is easy: {a„} are 
the moments of some measure supported on R if and only if all the Hankel 
determinants |ai+j|™ _Q, m = 0, 1,..., are positive. The unicity (usually called 
determinacy) depends on the behavior of the orthogonal polynomials (2.2) de¬ 
fined from the moments = ax+j . In fact, there are different measures with 
the same moments aj if and only if there is a non-real zq with |p„(zo) P < oo, 
which in turn is equivalent to |p„(z)p < oo for all z G C. Furthermore, the 
Cauchy transforms of all solutions n of the moment problem have the parametric 
form 

C di'{x) C{z)F{z) -\- A{z) 

J z-x ^ D{z)F{z) + B{z)' 

where F is an arbitrary analytic function (the parameter) mapping the upper 
half plane C+ into C+ U {oo}, and A, B, C and D have explicit representations 
in terms of the first and second kind orthogonal polynomials and Qn '■ 

Mz) = zJ2n^ri{0)qniz); B{z) = -1 z (0)p„(z); 

C'(z) = l-hzX;„Pn(0)g„(2); D{z) = zJ2nPni^)Pniz). 

For all these results see Akhiezer [2], and for an operator theoretic approach 
to the moment problem see Simon [78] (in particular, Theorems 3 and 4.14). 

Jacobi matrices and spectral theory of self-adjoint operators 

Tridiagonal, so called Jacobi matrices 


/ bo 

Oo 

0 

0 • 

•• \ 

ao 

bi 

ai 

0 • 


0 

ai 

b2 

02 • 


0 

0 

02 

62 


V : 




•• / 
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with bounded a„ > 0 and bounded real define a bounded self-adjoint operator 
J in I 2 , a so called Jacobi operator. These are the discrete analogues of second 
order linear differential operators of Schrodinger type on the half line. Every 
bounded self adjoint operator with a cyclic vector is a Jacobi operator in an 
appropriate base. 

The formal eigen-equation Jtt = Att is equivalent to the three-term recur¬ 
rence 

O^n — l'^n—l “t” ^n'^n “t” CLn'^n+1 — AtT^, Tl = 1,2,..., 

boTTg J- OoTTl = AtTo, TTq = 1. 


Thus, 7r„(A) is of degree n in A. 

By the spectral theorem, J, as a self-adjoint operator having a cyclic vector 
((1, 0,0,...)), is unitarily equivalent to multiplication by x in some with 

some probability measure /i having compact support on the real line. This fj, 
is called the spectral measure associated with J (and with its spectrum). More 
precisely, if Pn{x) = p„(/x,a:) are the orthonormal polynomials with respect 
to /X, and U maps the unit vector e„ = (0,..., 0,1, 0,...) into then U 
can be extended to a unitary operator from I 2 onto T^(/x), and if Sf{x) = 
xf{x) is the multiplication operator by x in Lp'^p), then J = U~^SU. The 
recurrence coefficients for p„(/i, x) are precisely the a„’s and 6„’s from the Jacobi 
matrix, i.e., Pn(x) = C7r„(x) with some fixed constant c. These show that Jacobi 
operators are equivalent to multiplication by x in (p) spaces if the particular 
basis {pn(M,-)} is used (see e.g. Deift [20, Ch. 2]). 

The truncated n x n matrix 


k bo ao 0 

oo bi ai 

0 oi 62 


0 

0 

02 




\ 0 0 


0 ^n—2 ^n—1 / 


has n real and distinct eigenvalues, which turn out to be the zeros of p„, i.e., 
the monic polynomial Pn{z)/Hn is the characteristic polynomial of J„. 


Quadrature 

For a measure /x, an n-point quadrature (rule) is a sequence of n points 
Xn,i, ■ ■ ■ ,Xn,n and a sequence of associated numbers Xn,i, ■ ■ ■ , ^n,n- It is ex¬ 
pected that 

/ n 

fdp ~ A„,fc/(Xn,fc) 
fc=l 

in some sense for as large a class of functions as possible. Often the accuracy of 
the quadrature is measured by its exactness, which is defined as the largest m 



Orthogonal Polynomials 


87 


such that the quadrature is exact for all polynomials of degree at most m, i.e., 
m is the largest number with the property that 


\r, 

k=l 


,kK 


for all 0 < j < m. 


For y, with support on the real line and for quadrature based on n points this 
exactness m cannot be larger than 2n — 1, and this optimal value 2n — 1 is 
attained if and only if Xn,i, • ■ •, Xn,n are the zeros of the orthonormal polynomial 
Pn{tk,x) corresponding to /i, and the so called Cotes numbers \n,k are chosen 
to be 


^n^k — ^n,k) — 


j=0 


Pj{pL,Xn,k)'‘ 


where Xn{p,z) is the Christoffel function (3.2) associated with p. 
See Szegd [91, Ch. XV]. 


Random matrices 


Some statistical-mechanical models in quantum systems use random matrices. 
Let Tin be the set of all n x n Hermitian matrices M = and let there 

be given a probability distribution on 7t„ of the form 

Pn{M)dM = D-^ exp(-nTr{V(M)})(iM, 

where fo(A), A € R, is a real-valued function that increases sufficiently fast at 
infinity (typically an even polynomial in quantum field theory applications), 
Tr{i7} denotes the trace of the matrix i7, 

n 

dM — dnik^k diRmkj d'^mkj 

k—1 k<j 


is the “Lebesgue” measure for Hermitian matrices, and Dn is a normalizing 
constant so that the total integral of Pn{M)dM is one. 

Every matrix M G has n real eigenvalues which carry physical informa¬ 
tion on the system when it is in the state described by M. The quantity 

^ ^ #{eigenvalues in D} 

n 

is the random variable that equals the normalized number of eigenvalues in the 
interval D. This model is known as the unitary ensemble associated with V. 

Let pj{vP,x) be the orthonormal polynomials with respect to the varying 
weight w'^{x), w{x) = exp(—V(x)). Then the joint probability density of the 
eigenvalues can be written in the form 


p*_i(w", Aj)w"/^(Aj) 
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where (i„ is a normalizing constant built up from the leading coefficients of the 
With the so called weighted reproducing kernel 

n — 1 

Knit,s) = 

0=0 

it can also be written in the form 

In particular, for the expected number of eigenvalues in an interval D we have 

EN„{B) = [ dX, 

Jd n 

where 1/Kn{X, A) is known in the theory of orthogonal polynomials as the n-th 
(weighted) Christoffel function associated with the weight w", while the limit 
of the left hand side (as n ^ oo) is known as the density of states. 

See, e.g., Mehta [60] and Pastur-Figotin [68]. 


7 Some questions leading to classical orthogonal 
polynomials 

There are almost an infinite number of problems where classical orthogonal 
polynomials emerge. Let us just mention a few. 


Electrostatics 

Put at 1 and —1 two positive charges p and q, and with these fixed charges 
put n positive unit charges on ]—1,1] at the points a;i,...,a;„. On the plane 
the Coulomb force is proportional with the reciprocal of the distance, and so a 
charge generates a logarithmic potential field. Therefore, the mutual energy of 
all these charges is 

n 1 n 1 

.= "g 'g ITTiM ^ 

and the equilibrium problem asks for finding xi,..., Xn for which this energy is 
minimal. The unique minimum occurs (see Szegd [91, Section 6.7]) for the zeros 
of the Jacobi polynomial orthogonal with respect to the weight 

(1 — x)^^“^(l + 

There is a similar characterization of the zeros of Laguerre and Hermite 
polynomials, and even of more general orthogonal polynomials (for the latter 
see Ismail [39, Section 3.5]). 
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Polynomial solutions of eigenvalue problems 

Consider the eigenvalue problem 

(P d 

/(a^)^J/(a;) + gix)—y{x) + h{x)y{x) = Xy{x), 

where f,g,h are fixed polynomials and A is a free constant, and it is required that 
this have a polynomial solution of exact degree n for all n = 0, 1 ,.. for which 
the corresponding A and y{x) will be denoted by A„ and t/n(x), respectively. 
Bochner’s theorem from [16] states that, except for some trivial solutions of the 
form y{x) = aa;" + bx'^ and for some polynomials related to Bessel functions, 
the only solutions are (in all of them we can take h{x) = 0) 

• Jacobi polynomials (fix) = 1 — g{x) = (3 — a — x{a + /3 + 2), 

A„ = -n(n + a + /3 + 1)) 

• Laguerre polynomials {f{x) = x, g{x) = 1 + a — x, A„ = —n) and 

• Hermite polynomials Hn{x) (/(x) = 1, g{x) = —2x, A„ = —2n). 


Harmonic analysis on spheres and balls 

Harmonic analysis on spheres and balls in is based on spherical harmonics, 
i.e., harmonic homogeneous polynomials. In this theory, special Jacobi polyno¬ 
mials, the so called ultraspherical or Gegenbauer polynomials play a 
fundamental role - they are orthogonal with respect to the weight (1 — 

Let 5"*“^ be the unit sphere in and let be the restriction to of 
all harmonic polynomials Q(xi,..., x„) of d variables that are homogeneous of 
degree n, i.e., 

" ' g2 

= Q(Axi,...,Ax„) = A”(5(xi,...,x„), a > 0. 

k=i 

The dimension of is 

/n-l-d — lA /n-bd — 3A 

V d-1 j"V d-1 )’ 

and an orthogonal basis in it can be produced as follows. With p = x'^_^ + x^ 

let gs,o = p’^Pi°\xd-i/p) and = Xdp‘‘Ps^\xd-i/p). With rid = 0 ornd = 1 
consider all multiindices n = (ni, n 2 ,..., Ud) such that ni -b • • • -I- rid = n, and 
if for such a multiindex we define the function Hn(xi,..., Xd) as 

d-2 

n + • • • + xlr^Pi^^\xj{x^j -b • • • -b , 

i=i 
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then these constitute an orthogonal basis in (see e.g. Dunkl-Xu [24, p. 
35]). 

If X = (xi,..., x„) and (x, y) = XkVk is the inner product in R'^, then 
the reproducing kernel for these spherical polynomials is y)) in the 

sense that for all Q € TLn and for all x S S‘^~^ we have 

c„,d J Pi‘^~'^^^'^{{x,y))Qiy)da{y) = Q{x), 

where integration is with respect to surface area, and Cn^d is an explicit normal¬ 
izing constant (see e.g. Dunkl-Xu [24, p. 37]). 

As a result, Gegenbauer polynomials are all over the theory of spherical 
harmonics, as well as in the corresponding theory for the unit ball. 


Approximation theory 

In the literature, expansions of functions into classical orthogonal polynomial 
series are second only to trigonometric expansions, and numerous works have 
been devoted to their convergence and approximation properties, see e.g. Szego 
[91, Ch. XIII]. 

The Chebyshev polynomials cos(narccosx) are orthogonal on [—1,1] 
with respect to the weight w{x) = (1 — x^)“^/^. These directly correspond 
to trigonometric functions, and expansions into them have virtually the same 
properties as trigonometric Fourier expansions. But there are many other as¬ 
pects of approximation where Chebyshev polynomials appear. If one considers, 
for example, the best approximation on [—1,1] of x" in the uniform norm by 
polynomials P„_i(x) of smaller degree then the smallest error appears when 
x" — P„_i(x) = 2^“" cos(narccosx) is the monic n-th Chebyshev polynomial. 
Actually, monic Chebyshev polynomials minimize all U’{w), p > 0, norms 
among monic polynomials of a given degree. 

As we have seen in (3.1), the monic orthogonal polynomials Pn(/r)/K„ are 
the solutions to the extremal problem 

y |P„|^dp ^ min, (7.3) 

where the minimum is taken for all monic polynomials of degree n. This ex¬ 
tremal property makes orthogonal polynomials, in particular Chebyshev poly¬ 
nomials, indispensable tools in approximation theory. 

Lagrange interpolation and its various generalizations like Hermite-Fejer or 
Hermite interpolation etc. is mostly done on the zeros of some orthogonal poly¬ 
nomials. In fact, these nodes are often close to optimal in the sense that the 
Lebesgue constant increases at the optimal rate. In many cases interpolation on 
zeros of orthogonal polynomials has special properties due to explicitly calcula¬ 
ble expressions. Recall e.g. Fejer’s result that if P 2 n-i is the unique polynomial 
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of degree at most 2n — 1 that interpolates a continuous function / at the nodes 
of the n-th Chebyshev polynomial and that has zero derivative at each of these 
nodes, then P 2 n-i uniformly converges to / on [—1,1] as n —>■ oo. For the role 
of orthogonal polynomials in interpolation see the books Szabados-Vertesi [90] 
and Mastroianni-Milovanovic [55]. 


8 Heuristics 


In this section we do not state precise results. We just want to indicate some 
heuristics on the behavior of orthogonal polynomials. For the concepts below, 
as well as for a more precise form of some of the heuristics see the following 
sections, in particular Section 9. 

As we have seen, the monic orthogonal polynomials Pn(/r)/Kn minimize the 
L^(/i) norm in (7.3). Therefore, these polynomials try to be small where the 
measure is large, e.g. one expects the zeros to cluster at the support S{fi) of 
yL. The example of arc measure on the unit circle, for which the orthogonal 
polynomials are z", shows however, that this is not true (due to the fact that 
the complement of the support is not connected). The statement is true when 
the support lies on R or on some systems of arcs, and also in the general case 
when instead of the support one considers the polynomial convex hull of the 
support of yt (for the definition, see the next section): on any compact set 
outside the polynomial convex hull there can only be a fixed number of zeros of 
p„(/x) for every n. When the complement of S{yt) is connected and S{y) has no 
interior, then the distribution of the zeros shows a remarkable universality and 
indifference to the size of y. In many situations the distribution of the zeros 
is the equilibrium distribution of the support S{yi). When S{y,) = [—1,1], this 
means that under very weak assumptions the zero distribution is always the 
arcsine distribution dxjTT\/l — . 

The L'^iy) minimality of Pn{y)/Kn in the sense of (7.3) is something like 
minimality in the L°° norm on S{fj,). Therefore, Pn{y)/nn should behave like 
the monic polynomial T„ minimizing the L°° norm on S{yL) (so called Chebyshev 
polynomials for Since 

-loglr„(z)|= [ \og\z - t\diy„{t) 

n J 


where Vn has mass 1/n at each zero of T„, in the limit the behavior should be 
like 


U''{z) = / logjz - t\dn{t), 


( 8 . 1 ) 


where v is the probability measure on S{yi) for which the maximum of IP' on 
S{y) is as small as possible (this is the so called equilibrium measure of S{fj,)). 
More generally, if dv = dvn = w''^{x)dx is a varying weight in the specified way, 
then the same reasoning leads to a behavior like (8.1), but now ^ is a measure 
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for which the supremum of U''{z) + logw(z) is as small as possible (weighted 
equilibrium measure). 

Universal behavior can also be seen for the polynomials themselves. Usually 
they obey 

-\og\pn{^i,z)\-> gc\s(^l:){z,cG), z ^ S{p), (8.2) 

where gc\s(p)iz, oo) is the Green function with pole at infinity associated with 
the complement of the support. When the unbounded component of the com¬ 
plement of S'(/i) is simply connected, then in that component often there is a 
finer asymptotic behavior of Pn(/r) of the form 

Pn{z) ~ d„5^(z)$(z)”, z ^ S{p), (8.3) 

where $ is the mapping function that maps C \ S{p) conformally onto the 
outside of the unit disk, and is a function (might be called generalized Szegd 
function) that depends on p. Such a fine asymptotic is restricted to the simply 
connected case, see e.g. Section 10. 

Asymptotics of orthogonal polynomials have a hierarchy, and the different 
types of asymptotics usually require the measure to be sufficiently strong with 
different degree on its support. Consider first the case of compact support 
S{p). The weakest is n-th root asymptotics stating the behavior (8.2) for 
|p„(/i, z)|^/" outside the support of the measure. It is mostly equivalent to a 
corresponding distribution of the zeros, as well as asymptotical minimal behav¬ 
ior of It holds under very weak assumptions on the measure, roughly 

stating that the logarithmic capacity of the points where p' > 0 (derivative 
with respect to equilibrium measure), be the same as the capacity of S{p). Ra¬ 
tio asymptotics, i.e., asymptotic behavior oi pn+i{p, z)/pnip, z), is stronger, 
and is equivalent with asymptotics for the ratio of consecutive leading 

coefficients. It can only hold when C\S{p) (more precisely its unbounded com¬ 
ponent) is simply connected, and in this case it is enough that p' > 0 almost 
everywhere with respect to the equilibrium measure of the support of p (see 
Section 10). Finally, strong asymptotics of the form (8.3) needs roughly that 
log^' be integrable (Szego condition, see Section 10). 

All these are outside the support. On the support the orthogonal polyno¬ 
mials are of oscillatory behavior, and in the real case under smoothness as¬ 
sumptions on the measure often a so called Plancherel-Rotach type asymptotic 
formula 

Pn{p, x) ~ dng{x) sin{nh{x) -I- H(x)) 

holds, where (/, h, H are fixed functions. Here h{x) is directly linked with the 
zeros, h' /ti is precisely the distribution of the zeros. When S{p) = [—1,1] and 
the measure is smooth, then h{x) = arccoscc. 

When S{p) is not of compact support (like Laguerre, Hermite or Freud 
weights), then usually the zeros are spreading out, and one has to rescale them 
to [—1,1] (or to [0,1]) to get a distribution, which is mostly NOT the arcsine 
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distribution. In a similar fashion, various asymptotics hold for the polynomials 
only after a corresponding rescaling. 

The Christoffel function (3.2) or, what is the same, the square sum 


1 

An(/r,2;) 


\pkiti, z)\‘^, 

k=0 


behaves much more regularly than the orthogonal polynomials. Outside the 
support the behavior of A„ is what one gets from the heuristics above on the 
polynomials (after square summation). On S'(/r) the typical behavior of A„(/i, z) 
is like /i(D„(z)), where D„(z) is the disk about 2 with equilibrium measure 
1/n (equilibrium measure of the support S'(/r)). In particular, nA„(/i, z) tends 
pointwise to the Radon-Nikodym derivative of fj, with respect to the equilibrium 
measure. As a rule of thumb the estimate |p„(/i,a:)p < C/nXn{p,x) holds in 
many cases. 

If / G T^(m)> its Fourier expansion into {pn{p, •)} is 


fix) ~ ^ k-kPkix)', Cfc 

k=0 

The n-th partial sum has the closed form 


fpkdfi. 


f{t)Knix,t)dfi{t), 


Kn{x,t) = '^Pk{x)pk{t). 

k=0 


In the real case for the reproducing kernel K„(x,t) we have the Christoffel- 
Darboux formula 


d^n (Xj t) 


Kn Pn+lix)Pn{t) - Pnix)Pn+lit) 

X t 


which suggests a singular integral-type behavior for the partial sums. In general, 
Fourier expansions into orthogonal polynomials are sensitive to the weight (re¬ 
call e.g. Pollard’s theorem that Legendre expansions are bounded in 1,1] 
only for 4/3 < p < 4), but sometimes convergence properties are equivalent to 
those of a related trigonometric Fourier series (so called transplantation theo¬ 
rems, see e.g. Askey [11]). 


9 General orthogonal polynomials 

In this section p always has compact support S{p). For all the results below 
see Stahl-Totik [87] and the references therein. 
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Lower and upper bounds 

The energy V{K) of a compact set K is defined as the infimum of 


J log (9-1) 

where the infimum is taken for all positive Borel measures on K with total 
mass 1. The logarithmic capacity is then cap(itr) = For the leading 

coefficients of the orthonormal polynomials Pn(/r) we have 


1 

cap{S{^i)) 


< liminf 

n^oo 


(9.2) 


To get an upper bound we need the concept of carrier: a Borel set if is a 
carrier for /r if /i(C \ E) = 0. The capacity of a Borel set is the supremum of 
the capacities of its compact subsets, and the minimal carrier capacity 
associated with /i is the infimum of the capacities of all carriers. With this 

limsupKy”<—. (9.3) 

n— 


When cap(iF) is positive, then there is a unique measure v = lok minimizing 
the energy in (9.1), and this measure is called the equilibrium measure of 
K. Green’s function gc\K{z,oo) with pole at infinity of C \ iF can then be 
defined as 

SO\k(..oo) = (f-'l) 

We have for all g (with cap(S'(^)) > 0) the estimate 

liminf-log |p„(/x, z) I > gc\s{fi){z,oo) (9.5) 

n—^oo 77, ' ^ ^ 


locally uniformly outside the convex hull of S{fj,), while in the convex hull but 
outside the so called polynomial convex hull Pc(5'(/x)) (for the definition see 
below) (9.5) is true quasi-everywhere (i.e., with the exception of a set of zero 
capacity). The same is true on the outer boundary of S{g,), which is defined 
as the boundary dfl of the unbounded component Q of the complement C\S{p,), 
namely for quasi-every z G dfl 

liminf |p„(^,z)I> 1. 

n —>oo 


The minimal carrier Green function gfi{z, oo) is the supremum for all carri¬ 
ers E of the Green function of G \ if, where the latter is defined as the infimum 
of gc\K for all compact subsets K of E. With this, 

limsup-log|p„(/r, z)|fo” < g^(z, oo) 

n—^oo ^ 


(9.6) 
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locally uniformly on the whole plane. 

All these estimates are sharp. 

When the bounds in (9.2) and (9.3) coincide we have convergence for 
and these bounds coincide precisely when the bounds in (9.5) and (9.6) do so. 


Zeros 


The zeros of p„(/r) always lie in the convex hull of the support of the 

measure y.. This is a consequence of the extremal property (3.1) of orthogonal 
polynomials. In fact, if there was a zero zq of p„(/r, z) outside the convex hull 
of the support, then we could move that zero towards the convex hull (along 
a line segment that is perpendicular to a line separating zq and During 

this move the absolute value of the polynomial decreases at all points of S'(/r) 
and hence so does its T^(/r) norm, but that is impossible by (3.1). 

To say somewhat more on the location of zeros we need the concept of the 
polynomial convex hull. When Q is the unbounded component of the com¬ 
plement C \ S{y,), then Pc(5'(/i)) = C \ 11 is called the polynomial convex 
hnll of S'(/i) (it is the union of S'(/i) with all the “holes” in it, i.e., with the 
bounded components of C \ S{fi)). Now the zeros cluster on Pc(5'(/r)) in the 
sense that for any compact subset AT of O there is a number Nk independent 
of n, such that Pn(/r) can have at most Nk zeros in K. The proof of this is 
based on the following lemma: Let V, S C C be two eompaet sets. If V and 
Pc(5') are disjoint, then there exist a < 1 and m G N such that for arbitrary 
m points xi,, Xm G V there exist m points yi,... ,ym G C for which the 
rational function 


r-miz) := 

1=1 


z-Vj 

Z-Xj 


(9.7) 


has on S a sup-norm satisfying 


||r^||s<a. (9.8) 

Taking this for granted, assume that is a compact set contained in H. We 
apply the lemma with S = S{y,), and let a < 1 and m G N be the numbers 
in the lemma. Let us assume that Pnit^', z) has at least m zeros xi,..., Xm on 
V. By the lemma there exist m points yi,... ,ym G C such that the rational 
function defined as in (9.7) by the points xi,..., Xm and yi,. .. ,ym satisfies 
the inequality (9.8). With r^ we define the modified monic polynomial 

quiz) := rra{z)pn{pi.-, z) = z” -I- . . . , 

For the L^(/x) norm of this polynomial we have the estimate 
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which contradicts the minimality (3.1) of the monic orthogonal polynomial 
p„(/r; z). Hence, we have proved that z) has at most m — 1 zeros on V, as 
was stated. 

What we have said about the zeros can be sharpened for measures on the 
real line. For example, if p is supported on the real line, then Pc(S'(/r)) = S{p), 
and if iF is a closed interval disjoint from the support, then there is at most 
one zero in K. It was shown in Denison-Simon [23] that if xq G R is not in 
the support, then for some (5 > 0 and all n either or Pn+i has no zero in 
{xq — i5, xq + d). Note that if /r is a symmetric measure on [—1, —1/2] U [1/2,1], 
then p2n+r(0) = 0 for all n, so the result is sharp. 

Any isolated point in the support that lies on the outer boundary attracts 
precisely one zero. Let zq be an isolated point of S'(/r), such that its distance 
from the convex hull of S{p) \ {zq} is i5 > 0. Then has at most one zero in 
the disk {jz — zq] < <5/3} (Simon [80, Section 8.1]). It is also clear that for any 
symmetric measure p with S{p) = [—1, —1/2] U {0} U [1/2,1] the polynomials 
P2n{p) have 2 zeros near 0, so the result is sharp (in this case 5 = 0). Moreover, 
if p lies on the unit circle, then there exist two positive constants C and a and 
a zero Zn of such that \zn — zq\ < Ce““". 

In general, each component of the polynomial convex hull consisting of more 
than one point attracts infinitely many zeros: if 7 is a Jordan curve in H such 
that S{p) n 7 is infinite, then the number of zeros of Pn that lie inside 7 tends 
to infinity (Saff-Totik [76]). Mass points of p do not necessarily attract zeros 
(above we have mentioned that they do if they lie on the outer boundary). 
In fact, it was shown in Saff-Totik [76] that if p is the measure on the unit 
circle given by the density function sin^( 0 / 2 ), then for any measure cr that is 
supported in the open unit disk there is a A > 0 such that all zeros of the n-th 
orthogonal polynomials with respect to p = p + Xa tend to the unit circle as 
n 00 . 

Next put a unit mass at every zero of Pnip) (counting multiplicity). This 
gives the so called counting measure on the zero set. Zero distri¬ 
bution amounts to finding the limit behavior of The normalized arc 

measure on the unit circle (for which pn{p, z) = z") shows that if the interior 
of the polynomial convex hull Pc(S'(/i)) is not empty, then the zeros may be 
far away from the outer boundary dfl, where the equilibrium measure x!s{fi) is 
supported. Thus, assume that Pc(S'(/i)) has empty interior and also that there 
is no Borel set of capacity zero and full /i-measure, i.e., the minimal carrier 
capacity is positive (the = 0 case is rather pathological, almost anything 
can happen with the zeros then). In this case 



cap(S'(^)) 


(9.9) 


lim-z.p„(^) =u;s(^) 


if and only if 
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in weak* sense, i.e., asymptotically minimal behavior of (see (9.2)) is equiv¬ 
alent to the fact that the zero distribution is the equilibrium distribution. In a 
similar way, asymptotic maximal behavior (see (9.3)), i.e., 

lim 4 /” = — (9.10) 

n^oo 


holds precisely when 

lim = uJf,, 

where ujfj, is the so called minimal carrier equilibrium measure, for which a 
representation like (9.4) is true, but for the minimal carrier Green function g^. 

Regularity 

(9.9) is called regular limit behavior, and in this case we write /i G Reg. 
Thus, the important class Reg is defined by the property (9.9). y. G Reg is 
equivalent to either of 

• lim„^oo bn(M. 2 )y” = exp( 5 c\s(p)( 2 ,oo)), z ^ Con{S{y,)) 

• limsup„^oo bn(M, .z)^" = 1 for quasi-every z; G dfl. 

If n is a regular set with respect to the Dirichlet problem, then /i G Reg is 
equivalent to either of 

• lim„^^ lbn(M)ll™”s(^) = 1 

• For any sequence {Pn} of polynomials of degree n = 1,2,... 

.. / H-PK||sup,5(At) 

WPnhH,) 

The last statement expresses the fact that in the n-th root sense the (/i) and 
L°° norms (on S'(/i)) are asymptotically the same. 

All equivalent formulations oi n G Reg point to a certain “thickness” of 
/i on its support. Regularity is an important property, and it is desirable to 
know “thickness” conditions under which it holds. Several regularity criteria 
are known, e.g. either of the conditions 

• all Borel sets B C S{fj,) with full measure (i.e with = fj,{S{y,))) have 

capacity cap(i3) = cap(5'(/x)), i.e., = cap(5'(/x)) or 

• dfj,/dujs{fj,) > 0 (Radon-Nikodym derivative) a;s(p)-almost everywhere 
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is sufficient for /i G Reg. Regularity holds under fairly weak assumptions on 
the measure, e.g. if 5'(/x) = [0,1], and 

lim inf r log /i([cc — r, a; + r]) >0 

for almost every x G [0,1] (i.e., if /i is not exponentially small around almost 
every point), then ^ G Reg. 

No necessary and sufficient condition for regularity in terms of the size of 
the measure /i is known. The only existing necessary condition is for the case 
S'(^) = [ 0 , 1 ], and it reads that for every ry > 0 

lim cap ({x | /r([x — 1 /n, x + l/n]) > ^ 

(here 1/4 is the capacity of [0,1]). 


10 Strong, ratio and weak asymptotics 


Strong asymptotics 

Let /i be supported on [—1,1] and suppose that the so called Szego condition 


1-1 Vl - 


dt > —oo 


( 10 . 1 ) 


holds, where fj! is the Radon-Nikodym derivative of /x with respect to linear 
Lebesgue measure. Note that this condition means that the integral is finite, 
for it cannot be oo. It expresses a certain denseness of /i, and under this con¬ 
dition G. Szego proved several asymptotics for the corresponding orthonormal 
polynomials p„(/x). This theory was developed on the unit circle and then was 
translated into the real fine. The Szego function associated with /i is 


:= exp 




\ogn'{t) dt ^ 
Z - t y/l J 


( 10 . 2 ) 


and it is the outer function in the Hardy space on C \ [—1,1] with boundary 
values \Dfj,{x)\‘^ = fi'{x). Outside [—1,1] the asymptotic formula 

p„(/x, z) = (1 -p o(l))—^(z -k \/z2 - l)”T)^(z)"^ (10.3) 

V ZTT 

holds locally uniformly. In particular, the leading coefficient of Pni^J-) is of 
the form 


Kn = il+ 0(1))—= exp 


/-I 


log 


J_, 


dt 


(10.4) 
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If c?/i(a;) = w{x)dx and h(t) = w(cost)sint satisfies a Dini-Lipshitz condition 

+ £> 0 , 

then with 


1 r logw(0 - logw(x) 


we have uniformly on [—1,1] 

(1 — ^ cos ^(n + i) arccosx + ru,(x) — 

+ 0((logn)“‘'). 


For all these results see Szegd [91], Chapter 6. The Szego condition is also 
necessary for these results, e.g. an asymptotic formula like (10.3) and (10.4) is 
equivalent to (10.1). 


Ratio asymptotics 

If one assumes weaker conditions then necessarily weaker results will follow. A 
large and important class of measures is the Nevai class M{b,a) (see Nevai 
[62]), for which the coefficients in the three-term recurrence 

XPn{x) = anPn+l{x) + bnPn{x) P an-lPn-l{x) 

satisfy a„ ^ a, bn ^ b. This is equivalent to ratio asymptotics 

p„+i(^:) z-b+ yj{z-by - 
lim -= --- 

n^oo Pn{z) 2 

for large z (actually, away from the support of p), and the monograph Nevai 
[62] contains a very detailed treatment of orthogonal polynomials in this class. 
It is also true that if the limit of Pn+i{z) / Pn(z) exists at a single non-real z, 
then p G M{b, a) for some a, b (Simon [79]). 

The classes M{b,a) are scaled versions of each other, and the most im¬ 
portant condition ensuring M(0,1/2) is given in Rakhmanov’s theorem from 
[75]: if p is supported in [—1,1] and p' > 0 almost everywhere on [—1,1], 
then p G M(0,l/2). Conversely, Blumenthal’s theorem from [15] states that 
p e M(0,1/2) implies that the support of p is [—1,1] plus at most count¬ 
ably many points that converge to ±1. Thus, in this respect the extension of 
Rakhmanov’s theorem given in [22] by Denisov is of importance: ii p' > 0 almost 
everywhere on [—1,1] and outside [—1,1] the measure p has at most countably 
many mass points converging to ±1, then p G M{0, 1/2). However, M(0,1/2) 
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contains many other measures not just those that are in these theorems, e.g. in 
Delyon-Simon-Souillard [21] a continuous singular measure in the Nevai class 
was exhibited, and the result in Totik [92] shows that the Nevai class contains 
practically all types of measures allowed by Blumenthal’s theorem. 

Weak and relative asymptotics 

Under Rahmanov’s condition supp(/r) = [—1,1], /r' > 0 a.e., some parts of 
Szegd’s theory can be proven in a weaker form (see e.g. Mate-Nevai-Totik 
[57, 58]). In these the Turan determinants 

r„(x) :=pf^(x) - Pn-lix)pn+lix) 

play a significant role. In fact, then given any interval D C (—1,1) the Turan 
determinant T„ is positive on D for all large n, and Tn(x)~^dx converges in the 
weak* sense to dp on D. Furthermore, the absolutely continuous part p' can 
be also separately recovered from Tn- 

lim [ Tn{x)p'{x) -(1 — dx = 0. 

n—yoc J TT 

Under Rahmanov’s condition we also have weak convergence, for example, 

lim [ f{x)pl{x)p'{x)dx = - [ —==dx (10.5) 

J TT J-i Vl - x^ 

for any continuous function /. Pointwise we only know a highly oscillatory 
behavior: for almost all x € [—1,1] 

limsupp„(x) > -(^'(a;))“^/^(l - 

n^oo TT 

liminfp„(x) < — — {p'{x))~^^^{l — 

n—TT 

and if En{e) is the set of points x G [—1,1] where 

\Pn{.x)\ > (1 + £)-{p'{x))~^/'^{l - 
TT 

then \En{£)\ —> 0 for all e > 0. However, it is not true that the sequence 
{pn{p, a^)} is pointwise bounded, since for every e > 0 there is a weight function 
w > 1 on [—1,1] such that Pn{Q)lrd-^‘^~’^ is unbounded (see Rakhmanov [74]). 

Simon [79] extended (10.5) by showing that if the recurrence coefficients 
satisfy b„ b, a 2 n+i a' and a 2 n o", then there is an explicitly calculated 
measure p depending only on 6, o', a" such that 

lim f f{x)pl{x)p'{x)dx = [ f{x)dp{x) (10.6) 
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for any continuous function /, and conversely, if (10.6) exists for f(x) = x, x'^, x^, 
then bn b, a 2 n+i cl' and a2„ ^ a" with some b, a', a”. 

For measures in Nevai’s class, part of Szego’s theory can be extended to 
relative asymptotics, i.e., when sequences of orthogonal polynomials corre¬ 
sponding to two measures are compared. Here is a sample theorem: let a be 
supported in [—1,1] and in Nevai’s class M(0,1/2), and let df3 = gda, where g 
is a function such that for some polynomial R both Rg and R/g are Riemann 
integrable. Then 


lim 

n—*oo 


Pn{P,z) 

Pn{a,z) 




uniformly on C away from [—1,1], where Dg is Szegd’s function with respect to 
the measure g{x)dx. 


Widom’s theory 

Szegd’s theory can be extended to measures lying on a single Jordan curve 
or arc J (see Kaliaguine [40] where also additional outside lying mass points 
are allowed), in which case the role of z + \/z^ — 1 in (10.3) is played by the 
conformal map 4) of C \ J onto the exterior of the unit disk, and the role of 
2" in (10.4) is played by the reciprocal of the logarithmic capacity of J (see 
Section 9). Things change considerably if the measure is supported on a set J 
consisting of two or more smooth curve or arc components Ji,, Jm. A general 
feature of this case is that «:„cap(J)” does not have a limit, its limit points fill 
a whole interval (though if some associated harmonic measures are all rational 
then the limit points may form a finite set). The polynomials themselves have 
asymptotic form 

^ = cap(J)"$(0)"(F„(0) + o(l)) 

Kn 

uniformly away from J, where $ is the (multi-valued) complex Green function of 
the complement C \ J, and where Fn is the solution of an L^-extremal problem 
involving analytic functions belonging to some class r„. The functions F in 
Tji are determined by an FI^ condition plus an argument condition, namely 
if the change of the argument of <I> as we go around Jk is 7fe27r modulo 27r, 
then in r„ we consider functions whose change of the argument around Jk is 
—n7fe27r modulo 27r. Now the point is that these function classes r„ change 
with n, and hence so does Fn, and this is the reason that a single asymptotic 
formula like (10.4) or (10.3) does not hold. The fundamentals of the theory 
were laid out in H. Widom’s paper [97]; and since then many results have been 
obtained by F. Peherstorfer and his collaborators, as well as A. I. Aptekarev, J. 
Geronimo, S. P. Suetin and W. Van Assche. The theory has deep connections 
with function theory, the theory of Abelian integrals and the theory of elliptic 
functions. We refer the reader to the papers Aptekarev [7], Geronimo~Van 
Assche [31], Peherstorfer [69]-[72] and Suetin [88]-[89]. 
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Asymptotics for Christoffel functions 

The Christoffel functions 

n 

k=0 

behave somewhat more regularly than the orthogonal polynomials. In Mate- 
Nevai-Totik [59] it was shown that if p is supported on [—1,1], it belongs to 
the Reg class there (see Section 9) and log/i' is integrable over an interval 
I C [—1,1], then for almost all x G I 

lim nXn{p,x) = tta/ 1 — x^p'{x). 

n—yoo 


This result is true (see Totik [93]) in the form 

lim nXn{p,x) = -- ^ , a.e. x G / 

«Wsupp(^) (X) 

when the support is a general compact subset of R, p G Reg and logp' G L^{I). 

Often only a rough estimate is needed for Christoffel functions, and such a 
one is provided in Mastroianni-Totik [56]: if p is supported on [—1,1] and it is 
a doubling measure, i.e., 

p{21) < Lp{I) 

for all / C [—1,1], where 21 is the twice enlarged I, then uniformly on [—1,1] 




An(x) 


x/l — 1 

n 


11 Recurrence coefficients and spectral measures 

Let p he a unit measure of compact support on the real line, and 

XPn{x) = anPn+l{x) + bnPn{x) + a„_ip„_l(x), 

the recurrence relation for the corresponding orthogonal polynomials. We have 
already mentioned in Section 6 that p is the spectral measure for the Jacobi 


matrix 
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/ 


and this gives a one-to-one correspondence between unit measures with compact 
(and infinite) support on the real line and Jacobi operators with bounded en¬ 
tries. Every such Jacobi operator is a bounded self-adjoint operator on I 2 , hence 
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operator theory and orthogonal polynomials meet at this point, and techniques 
and questions from both areas are relevant. Information between the measure 
H and the recurrence coefficients might be called the spectral analysis of or¬ 
thogonal polynomials. In this, a special role is played by the Chebyshev case 
a„ = 1/2 and = 0, where the Jacobi matrix is denoted by Jq. If the sequences 
{a„} and {bn} have limits then we may assume a„ ^ 1/2, ^ 0 (this is just 

rescaling, and this was the Nevai class M(0,1/2)), hence in this case the Jacobi 
operator is a compact perturbation of Jq, and one of the main questions of the 
theory is how properties of J — Jq are reflected in the spectral measure 

We have already mentioned in Section 10 that /r G M(0,1/2) implies that 
the support of ^ is [—1,1] plus some additional mass points converging to ±1 
(called Blumenthal’s theorem in orthogonal polynomials; it is a special case of 
Weyl’s theorem in operator theory on the invariance of the essential spectrum 
under compact perturbation). Conversely, if the support of /i is [—1,1] plus 
some additional mass points converging to ±1 and /i'(x) > 0 for almost all 
X G [—1,1], then /X G M{a,b) (Denisov [22], Rakhmanov [75], Nevai-Totik [63]). 
No spectral characterization of n G M(0,1/2) is known; this important class 
seems to contain all sorts of measures. For example, if v is any measure with 
support [—1,1] then there is a yx G M(0,1/2) which is absolutely continuous with 
respect to v. In particular, M(0,1/2) contains discrete measures, continuously 
singular measures or measures that are given by a continuous density which is 
positive on a set of measure < e. 

Strengthening the condition G M(0,l/2) can be done in several ways. 
After numerous works in the subject by Szegd, Shohat, Geronimus, Krein, Kol¬ 
mogorov and others, a complete characterization for J — Jq being a Hilbert- 
Schmidt operator was given in [43] by R. Killip and B. Simon (note that fj, is 
assumed to have total mass 1): 

^(a„ - 1/2)2+ ^62 ^ oo 

n n 

if and only if the following conditions hold: 

(i) the support of /i is [—1,1 ] plus some additional mass points converging 

to ±1, 

(ii) if /x' is the absolutely continuous part of ^ on [—1, 1], then 



(iii) for the mass points lying outside [—1,1] we have 

^|i?+-l|3/2+^|£:- + l|3/2<oo. 

J J 
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It was also shown in Killip-Simon [43] that if J — Jq is trace class, i.e., 


\an - 1 / 2 | + < oo, 


then Szego’s condition 



log ft) 
VT^ 


dt > —oo 


( 11 . 2 ) 


holds. The conclusion is also true if /i G M(0,1/2), for the mass points 
lying outside [— 1 , 1 ] we have 




and 

limsup(2"ai • • • a„) > 0. (11-3) 

n 

If the support of /i is contained in [—1,1], then Szegd’s condition automatically 
holds if (11.3) is true. Actually, when supp(/x) = [—1,1] then Szegd’s condition 
( 11 . 2 ) is equivalent to ( 11 . 1 ) and to the (conditional) convergence of the series 
- 1 / 2 ) and 

There is also an extended theory of orthogonal polynomials with several dif¬ 
ferent applications when the recurrence coefficients do not converge, but they are 
asymptotically periodic in the sense that for some k all the sequences {akn+j)'?!^:^i 
and {bkn+j)^=i, j = 1, • ■ •, fc converge. These are related to so called sieved or¬ 
thogonal polynomials and to orthogonal polynomials generated by polynomial 
mappings. In this case the essential support of the spectral measure lies on 
several intervals. There are numerous papers on this subject by M. E. H. Is¬ 
mail, N. A. Al-Salam, J. A. Gharris, J. Wimp, J. Bustoz, J. Geronimo, W. Van 
Assche, F. Peherstorfer, R. Steinbauer, N. I. Akhiezer, B. P. Osilenker and oth¬ 
ers; see e.g. Gharris-Ismail [17], Geronimo-Van Assche [31], Peherstorfer [70], 
Peherstorfer-Steinbauer [73], Akhiezer [1] for details and for further references. 

12 Exponential and Freud weights 

These are weight functions of the form where x is on the real line or 

on some subinterval thereof. For simplicity we shall first assume that Q is even. 
We get Preud weights when Q{x) = jxj“, a > 0, a; G R, and Erdos weights 
if Q tends to infinity faster than any polynomial as ja;j ^ oo. G. Freud started 
to investigate these weights in the sixties and seventies, but they independently 
appeared also in the Russian literature and in statistical physics. One can 
safely say that some of Freud’s problems and the work of P. Neva! and E. A. 
Rahmanov were the primary cause of the sudden revitalization of the theory of 
orthogonal polynomials since the early 1980’s. In the last 20 years D. Lubinsky 
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with coauthors have conducted systematic studies on exponential weights, see 
e.g. Levin-Lubinsky [46, 45], Lubinsky [49], Lubinsky-Saff [50], Van Assche 
[94]; we should mention the names E. Levin, E. B. Saff, W. Van Assche, E. A. 
Rahmanov and H. N. Mhaskar. In the mid 1990’s a new stimulus came from 
the Riemann-Hilbert approach that was used together with the steepest descent 
method by P. Deift, T. Kriecherbauer, K. T.-R. McLaughlin, S. Venakides and 
X. Zhou ([20]) to give complete asymptotics when Q is analytic. 

One can roughly say that because of the fast vanishing of the weight around 
infinity, things happen on a finite subinterval [—a„,a„] (depending on the de¬ 
gree of the polynomials), and on [—a„, a„] techniques developed for [—1,1] are 
applied. For Freud weights one can also make the substitution x and 

go to orthogonality with respect to the varying weight , in which case 

things are automatically reduced to a finite interval which is the support of a 
weighted energy problem. 

The On are the so called Mhaskar-Rahmanov-Saflf numbers defined by 


2 1'^ a„tQ'{a„t) 
TT Jo Vl - 


dt. 


( 12 . 1 ) 


The zeros oipn{w^), w(x) = exp(—Q(x)) are spreading out and the largest zero 
is very close to a„, which tends to oo. 

To describe the distribution of the zeros and the behavior of the polynomials 
one has to make appropriate contractions. Let us consider first the case of Freud 
weight iu(x) = exp(—ja;j“), and let pn be the n-th orthogonal polynomial with 
respect to (on (—00,00)). In this case 







Thus, for the largest zero Xn,n we have ^ 7 q as n ^ 00, and to 

describe zero distribution we divide (contract) all zeros x„^i by These 

contracted zeros asymptotically have the Ullman distribution 


dpw{t) 


a 


,a — 1 


dt TT 

This measure pw minimizes the weighted energy 


du, t G [—1,1]. 


log 


1 


X — t\ 


+ 2 / Qdfi 


( 12 . 2 ) 


(12.3) 


among all probability measures compactly supported on R. It is a general 
feature of exponential weights that the behavior of zeros of the polynomials is 
governed by the solution of a weighted energy problem (weighted equilibrium 
measures, see Saff-Totik [77]). If is the leading coefficient of i.e., Pn{z) = 
KnZ'^ + ■ ■ ■, then (Lubinsky-Saff [50]) 

lim K„7ri/22-"e-"/“n("+i/2)/a ^ 
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and we have 


lim b„(ni/“7a2)r/” 

n^oo 

= exp ^log \z + \/z2 _ i| + Re^ du^ 

locally uniformly outside [—1,1]. This latter is so called n-th root asymptotics, 
while the former is strong asymptotics. Strong asymptotics for Pn{z) on different 
parts of the complex plane was given using the Riemann-Hilbert approach, see 
Deift [20] and Kriecherbauer-McLaughlin [44] and the references there. On the 
real line we have a Plancherel-Rotach type formula 

exp(-n7“|a;|“)- 


v^l — a;- 


: COS 


1 -r TT 

- arccosx + n7r/Xu,([a;, 1]) — — 


uniformly on any subinterval of (—1,1). 

Things become more complicated for non-Freud weights, but the correspond¬ 
ing results are of the same flavor. In this case the weight is not necessarily 
symmetric, but under some conditions (like Q being convex or xQ'{x) being in¬ 
creasing for a; > 0 and an analogous condition for a; < 0) the relevant weighted 
equilibrium measure’s support is an interval, and the definition of the Mhaskar- 
Rahmanov-Saff numbers a±„ is 


n 

0 



_ xQ'{x) 

^/{x- a-n){a„ 

_ Q'jx) 

^/ix- a-„){a„ 


x) 


dx, 


=dx. 

x) 


Now one solves the weighted equilibrium problem (12.3) for all measures p with 
total mass n, and if /i„ is the solution then [a_„,a„] is the support of and 
Hn/n will play the role of the measure pw from (12.2) above. 

The weight does not even have to be defined on all R, e.g. in [46] a theory 
was developed by Levin and Lubinsky that simultaneously includes far reaching 
generalizations of non-symmetric Freud, Erdos and Pollaczek weights such as 
(a) nonsymmetric Freud-type weights 


Q{x) 


|a;|“, x G [0, oo) 
|a;|^, x G (—oo, 0), 


(b) nonsymmetric Erdos weights such as 


Q{x) = 


exp/(|a;|“) - exp;(0), a;G[0,oo) 

expfc(|a:|^) - expfc(O), a:G(-oo,0) 
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with exp; is the ?-times iterated exponential function, or 

(c) nonsymmetric Pollaczek type weights that vanish fast at ±1 such as 


r exp;((l-a:) “)-exp;(l), a; G [0,1) 

1 expj,((l - - expj,(0), a;G(-l,0). 


In all cases the interval [a_„,a„] is where things happen, e.g. this is the 
shortest interval on which the supremum norm of a weighted polynomial is 
attained: 

llsup = ll'*^-fn||sup,[—o_„,a„] 

for all polynomials of degree at most n. These numbers a±„ are everywhere in 
the theory, e.g. 

sup \pnix)\w(x)\x - ~ 1. 


13 Sobolev orthogonality 

In Sobolev orthogonality we consider orthogonality with respect to an inner 
product 

j (13.1) 

where pk are given positive measures. There are several motivations for this 
kind of orthogonality. Perhaps the most natural one is smooth data fitting. The 
Spanish school around F. Marcellan, G. Lopez and A. Martinez-Finkelshtein 
has been particularly active in developing this area (see the surveys Marcellan- 
Alfaro-Rezola [51] and Martinez-Finkelshtein [53, 52] and the references therein). 

In this section let Qn{z) = z" -|- • • • denote the monic orthogonal polyno¬ 
mial with respect to the Sobolev inner product (13.1), and qnipk) the monic 
orthogonal polynomials with respect to the measure pk- 

Most arguments for the standard theory fail in this case, e.g. it is no longer 
true that the zeros lie in the convex hull of the support of the measures pk, 
fc = 0,1,..., r. It is not even known if the zeros are bounded if all the measures 
/ifc have compact support. Nonetheless, for the case r = 1, and G Reg 

(see Section 9) it was shown in Gautschi-Kuijlaars [29] that the asymptotic 
distribution of the zeros of the derivative Q'„ is the equilibrium measure coequEi , 
where Ei is the support of pi, i = 0,1 (which also have to be assumed to be 
regular). Furthermore, if, in addition, Eq C Ei, then the asymptotic zero 
distribution of Qn is ujeo- 

In general, both the algebraic and the asymptotic/analytic situation is quite 
complicated, and there are essentially two important cases which have been 
understood to a satisfactory degree. 

Case I: The discrete case. In this case /xq is some “strong” measure, e.g. from the 
Nevai class M{b, a) (see Section 10), and pi,... ,pk are finite discrete measures. 
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It turns out that then the situation is similar to adding these discrete measures 
to iiq (the new measure will also be in the same Nevai class), and considering 
standard orthogonality with respect to this new measure. For ecample, if r = 1, 
then 


n^oo + fil,z) 

holds uniformly on compact subsets of C \ supp(^o + Mi)- Thus, the Sobolev 
orthogonal polynomials differ from those of the measure but not more than 
what happens when adding mass points to ^o- 

In this discrete case the Q„’s satisfy a higher order recurrence relation, hence 
this case is also related to matrix orthogonality (see the end of the Section 16). 


Case II: The Szego case. Suppose now that are all supported on 

the same smooth curve or arc J, and they satisfy Szegd’s condition there (see 
Section 10). In this case the fc-th derivative of Qn satisfies, locally uniformly in 
the complement of J, the asymptotic formula 

lini = 1 

where $ is the conformal map that maps C \ J onto the complement of the 
unit disk. That is, in this case the measures /xq, ..., Mr-i do not appear in 
the asymptotic formula, only matters. The reason for this is the following: 
Q = Qn minimizes 

(0,Q) = E (13.2) 

among all monic polynomials of degree n, while q — qn-k{lkk) minimizes 

[ IgpdMfc 


among all monic polynomials of degree n — k. But the polynomial Q^n\t) = 
n{n — 1) • • • (n — fc + l)t"“^ + • • • is a monic polynomial times the factor n{n — 
1) • • • (n — fc + 1) ~ n^, and this factor is dominant for k = r, so everything else 
will be negligible. There are results for compensation of this factor which 
lead to Sobolev orthogonality with respect to varying measures. 

Under the much less restrictive assumption that /xq G Reg (see Section 9) 
and the other measures /Xfc are supported in the support E of /xq it is true 
(Lopez-Pijeira-Cabrera-Izquierdo [47]) that the asymptotic zero distribution of 
is the equilibrium measure uie for all k, 

lim ||(5i''^|lsup,is = cap(i^), 

n—*oci 
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and hence, away from the zeros in the unbounded component of the complement 
of if, we have 

lim |Q^(2)r/” = 

n—*00 

where gc\E is the Green function for this unbounded component. 

The techniques developed for exponential weights and for Sobolev orthog¬ 
onality were combined in Geronimo-Lubinsky-Marcellan [30] to prove strong 
asymptotics for Sobolev orthogonal polynomials when r = 1 and /iq = are 
exponential weights. 


14 Non-Hermitian orthogonality 

We refer to non-Hermitian orthogonality in either of these cases: 

• the measure /i is non-positive or even complex-valued and we consider 
with 

J Pn{z)z'^dfi = 0, k = 0,l,...,n-l, (14.1) 

• /r is again non-positive or complex-valued, or positive but lies on a complex 
curve or arc and orthogonality is considered without complex conjugation, 
i.e., 

J Pn{z)z’^dii = 0, k = 0,1,... ,n — 1. (14-2) 

More generally, one could consider non-positive inner products, but we shall 
restrict our attention to complex measures and orthogonality (14.2). 

As an example, consider the diagonal Fade approximant to the Gauchy 
transform 

/(.) = / fli) 

of a signed or complex-valued measure, i.e., consider polynomials and qn of 
degree at most n such that 

f{z)pn{z) - qn{z) = 0 ( 2 :"”"^) 

at infinity. Then satisfies the non-Hermitian orthogonality relation 

J Pn{x)x^dp{x) = 0, j = 0,1,... ,n-1. (14-3) 

In this non-Hermitian case even the Gram-Schmidt orthogonalization pro¬ 
cess may fail, and then is defined as the solution of the orthogonality condi¬ 
tion (14.1), resp. (14.2), which give a system of homogeneous equations for the 
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coefficients of p„. Thus, may have smaller degree than n, and things can get 
pretty wild with this kind of orthogonality. For example, in the simple case 

diJ,{x) = {x — C0S7rQ;i)(a; — C0S7rQ;2)(l — x'^)~^^^dx, x G [—1,1], 


with 0 < ai < 02 < 1 rationally independent algebraic numbers, the zeros 
of Pn from (14.3) are dense on the whole complex plane (compare this with 
the fact that for positive p all zeros lie in [— 1,1]). In Stahl [86] it was shown 
that it is possible to construct a complex measure /x on [—1,1], such that for an 
arbitrary prescribed asymptotic behavior some subsequence } will have this 
zero behavior. Nonetheless, the asymptotic distribution of the zeros is again the 
equilibrium distribution of the support of p under regularity conditions on p. 
For example, this is the case if 

• 1^1 belongs to the Reg class (see Section 9), and the argument of p, i.e., 
dfj,{t)/d\p\{t), is of bounded variation (Baratchart-Kiistner-Totik [12]), or 

• dp{x) = (/(arccosa;)(l — x^)~^^‘^dx, x G [—1,1], g is bounded away from 
zero and infinity, and satisfies \g{9 + 6) — g{0)\ < K\ logJ]”^"'^, or 

• /X is supported on finitely many intervals, the argument of p is uniformly 

continuous and for a(S) = inf2,gsupp(^) — S,x + i5]) the property 

lim^^o loga((5) = 0 holds (Stahl [86]). 

In [84]-[86] H. Stahl obtained asymptotics for non-Hermitian orthogonal poly¬ 
nomials even for varying measures and gave several applications of them to Fade 
approximation. When the measure g is of the form dg{x) = g{x){l — x^)~^^‘^dx, 
X G [—1,1], with an analytic g, for z G C \ [—1,1], a strong asymptotic formula 
of the form 


Pn{z) 

Kn 




exp 




VT^ 



(with the Szegd function (10.2)) was proved by J. Nuttall [66], [67], A. A. 
Gonchar and S. P. Suetin [38]. For a recent Riemann-Hilbert approach see 
the paper [10] by A. I. Aptekarev and W. Van Assche. A similar result holds 
on the support of the measure, as well as for the case of varying weights, see 
Aptekarev-Van Assche [10]. 


15 Multiple orthogonality 

Multiple orthogonality comes from simultaneous Fade approximation. It is a 
relatively new area where we have to mention the names of E. M. Nikishin, 
V. N. Sorokin, A. A. Gonchar and E. A. Rahmanov, A. I. Aptekarev, A. B. 
J. Kuijlaars, J. Geronimo and W. Van Assche (see the survey [95] by W. Van 
Assche and the references there and the paper Gonchar-Rakhmanov [37]). The 



Orthogonal Polynomials 


111 


analogues of many classical concepts and properties have been found, and also 
the analogues of the classical orthogonal polynomials are known, e.g. in the 
multiple Hermite case the measures are dfj,{x) = e~^ 

Asymptotic behavior of multiple orthogonal polynomials is not fully under¬ 
stood yet due to the interaction of the different measures. For the existing 
results see Aptekarev [8], Van Assche [95], Van Assche’s Chapter 23 in [39] and 
the references therein. 

Types and normality 

On R let there be given r measures ni,..., fir with finite moments and infinite 
support, and consider multiindices n = {ni,... ,nr) of nonnegative integers 
with norm |n| = ni + ■ ■ ■ + Ur- There are two types of multiple orthogonality 
corresponding to the appropriate Hermite-Pade approximation. 

In type I multiple orthogonality we are looking for polynomials Qnj of 
degree Uj — 1 for each j = 1 ,..., r, such that 

J x^Qnj{x)dfij{x) = 0, A: = 0,1 ,..., |n| - 2. 

j=i 

These orthogonality relations give |n| — 1 homogeneous linear equations for the 
|n| coefficients of the r polynomials Qnj, so there is a non-trivial solution. If the 
rank of the system is |n| — 1, then the solution is unique up to a multiplicative 
factor, in which case the index n is called normal. This happens precisely if 
each Qn,j is of exact degree Uj — 1. 

In type II multiple orthogonality we are looking for a single polynomial 
Pn of degree |n| such that 

J x’^Pn{x)dyLi{x) = 0, A: = 1 ,..., m — 1 

J X^Pn{x)dyLr{x) = 0, /c = 1, . . . , Ur — 1. 

These are |n| homogeneous linear equations for the |n| -I- 1 coefficients of Pn, 
and again if the solution is unique up to a multiplicative constant, then n is 
called normal. This is again equivalent to being of exact degree n. 

n is normal for type I orthogonality precisely when it is normal for type 
II, so we just speak of normality. This is the case, for example, if the iij’s 
are supported on intervals [aj,bj] that are disjoint except perhaps for their 
endpoints; in fact, in this case has Uj simple zeros on (aj,bj). Normality 
also holds if dfij = Wjdg. with a common fj, supported on some interval [a, 5], 
and for all mj < Uj, j = 1 ,... ,r, every non-trivial linear combination of the 
functions 

Wi{x),XWi{x), . . . , x'^'^~^Wi{x),W 2 {x),XW 2 {x), . . . , x'^^~^Wr{x) 
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has at most mi + • • • + m^ — 1 zeros on [o, b] (this means that these functions 
form a so called Chebyshev system there). In this case has \n\ — 1 zeros on 
[a, 6]. 


Recurrence formulae 


To describe recurrence formulae, let = (0,..., 1,..., 0) where the single 1 
entry is at position j. Under the normality assumption if P„ is the monic 
orthogonal polynomial, then for any k 


i=i 


Another recurrence formula is 

r 

— Pn+e^ (^) d” bji ^Pj^X^ + ^ ^ -— 

1=1 


where 7r(l),..., 7r(r) is an arbitrary, but fixed, permutation of 1, 2,..., r. The 
orthogonal polynomials with different indices are strongly related to one an¬ 
other, e.g. Pn+ej^{x) — Pn+e^{x) is a Constant multiple of P„(x). 

If dfXj = Wjdfx, then similar recurrence relations hold in case of type I or¬ 
thogonality for 

r 

Qnj (yX^Wj . 

1=1 

Also, type I and type II are related by a biorthogonality property: 


PnQmd^^ — 0 


except for the case when m = foi some k, and then the previous integral 

is not zero (under the normality condition). 

To describe an analogue of the Christoffel-Darboux formula let {nij} be 
a sequence of multiindices such that tjIq is the identically 0 multiindex, and 
rrij 1 1 coincides with except for one component which is 1 larger than the 
corresponding component of rUj. Set Pj = Pm^, Qj = Qm.^^ and with m = m„ 

:= J Pm{x)x’'—'>^diij{x), 

where {m)j denotes the j-th component of the multiindex m. Then (see Daems- 
Kuijlaars [19]), again with m = m„, 

«-i j^U) 

{x-y)'^ Pk{x)Qk{y) = Pm{x)Q^{y) - ^ —(^—Prn-e. {x)Q^+e. {y) ■ 

k—0 1 — 1 ^m — e. 

Thus, the left hand side depends only on m = m„ and not on the particular 
choice of the sequence mj leading to it. 



Orthogonal Polynomials 


113 


The Riemann Hilbert problem 

There is an approach (see Van Assche-Geronimo-Kuijlaars [96]) to both types 
of multiple orthogonality in terms of matrix-valued Riemann-Hilbert problem 
for (r -I- 1) X (r -I- 1) matrices Y = (Vy (^))[j= 0 ' 

If dfj,j{x) = Wjdx, then one requires that 

• y is analytic on C \ R, 

• if Y^{x) denote the limit oiY{z) as 2 ; ^ a; G R from the upper, respec¬ 
tively the lower, half plane, then we have Y~^(x) = Y~{x)S{x), where 

1 wi{x) W2{x) Wr(x) 

01 0 ••• 0 

00 1 ••• 0 

00 0 ••• 1 

• as z ^ 00 

■ zliil 0 0 

0 0 

r(,) = (^/ + oQ)) 0 0 Z-- 

0 0 0 

The first entry yii(z) is precisely the orthogonal polynomial of type II, and 
the other entries are also explicit in terms of the P„’s and Wj ’s (all other entries 
are either a constant multiple of Pn-e^. or a Cauchy transform of its multiple 
with Wj ). For type I orthogonality the transfer matrix is 

1 0 0 ••• 0 

—wi(a;) 1 0 ••• 0 

S{x) := -^" 2 ( 3 ;) 0 1 • • • 0 

— Wr{x) 0 0 • • • 1 

the behavior at infinity is of the form 

■ z-liil 0 0 • • • 0 

0 z”i 0 • • • 0 

y(,)= + 0 0 z- ... 0 

_ 0 0 0 ... z" 

and the multiple orthogonal polynomials Qn,j are Yijj.i/2'!Ti. 
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16 Matrix orthogonal polynomials 

In the last 20 years the fundamentals of matrix orthogonal polynomials have 
been developed mainly by A. Duran and his coauthors (see also the work [9] by 
A. I. Aptekarev and E. M. Nikishin). The theory shows many similarities with 
the scalar case, but there is an unexpected richness which is still to be explored. 

For all the results in this section see Lopez-Rodriguez-Duran [48] and Duran- 
Griinbaum [25] and the numerous references there. 


Matrix orthogonal polynomials 


An N X N matrix 


Pit) 


( Piiit) 


PlN{t) 


\ PNl{t) ■ • ■ PNNit) / 

with polynomial entries Pij{t) of degree at most n is called a matrix polyno¬ 
mial of degree at most n. Alternatively, one can write 

p{t) = + • • • + Co 


with numerical matrices Cn, ■ ■ ■ ,Cq of size N x N. 

The number t = a is called a zero of P if T’(a) is singular, and the mul¬ 
tiplicity of a is the multiplicity of a as a zero of detP(a). When the leading 
coefficient matrix C„ is non-singular, then P has nN zeros counting multi¬ 
plicity. 

From now on we fix the dimension to be N, but the degree n can be any 
natural number. / will denote the N x N unit matrix and 0 stands for all kinds 
of zeros (numerical or matrix). 

A matrix 

/ /Xii(t) ••• fiiNit) \ 


W{t) = 


\ PNlit) ■ • ■ PNNit) / 

of complex measures defined on (or part of) the real line is positive definite 
if for any Borel set E the numerical matrix W(E) is positive semidefinite. We 
assume that all moments of W are finite. With such a matrix we can define a 
matrix inner product on the space of N x N matrix polynomials via 


(p,g) = j p{t)dwit)Q*{t), 


and if (P, P) is nonsingular for any P with nonsingular leading coefficient, then 
just as in the scalar case one can generate a sequence {Pn}'^=o of matrix poly¬ 
nomials of degree n = 0,1,... which are orthonormal with respect to W: 


I p^{t)dw{t)p:^{t) 


0 if n ^ m 
I ii n = m, 
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and here P„ has nonsingular leading coefficient matrix. The sequence {P„} is 
determined only up to left multiplication by unitary matrices, i.e., if C/„ are 
unitary matrices, then the polynomials also form an orthonormal system 

with respect to W. 

Three-term recurrence and quadrature 

Just as in the scalar case, these orthogonal polynomials satisfy a three-term 
recurrence relation 

tPn{t) = An+lPn+l{t) + BnPn{t) + A^Pn-i{t), U > 0, (16.1) 

where An are nonsingular matrices, and Bn are Hermitian. Conversely, the 
analogue of Favard’s theorem is also true: if a sequence of matrix polynomials 
{Pn} of corresponding degree n = 0 , 1, 2,..., satisfy (16.1) with nonsingular An 
and Hermitian Bn, then there is a positive definite measure matrix W such that 
the Pn are orthonormal with respect to W. 

The three-term recurrence formula easily yields the Christoffel-Darboux for¬ 
mula: 

n — 1 

{w-z)Y, P:{z)Pk{w) = P:_,{z)AnPn{w) - P*n{z) A*nPn-l{w), 

from which for example it follows that 

P:_,{z)AnPn{z) - P:{z)A*nPn-l{z) = 0 , 

n — 1 

^ P:iz)Pk{z) = P:_i(z)H„P;(z) - P:(z)H:p;_i(z). 
k^O 

The orthogonal polynomials Qn of the second kind 

_ f Priit) - Pn{x) „_1 o 

Qn{i) — / , dW (x), Tl — 1,2,..., 

J t-X 

also satisfy the same recurrence and are orthogonal with respect to some other 
matrix measure. For them we have 

P*_,{t)AnQn{t) - P*{z)A*nQn-l{t) = I, 

and 

Q„(t)P:_l(t) - Pnm*n-lit) = A-K 

With the recurrence coefficient matrices An, Bn one can form the block 
Jacobi matrix 


/ Bo 

Ao 

0 

0 ••• \ 

O* 

Pi 

Al 

0 ••• 

o 

Al 

B2 

A 2 ■ ■ ■ 

o • • • 

o • • • 

Al 

P 2 ••• 
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The zeros of P„ are real and they are the eigenvalues (with the same multiplicity) 
of the A^-truncated block Jacobi matrix (which is of size nN). If a is a zero 
then its multiplicity p is at most N, the rank of Pn{o) is N — p, and the space 
of those vectors v for which Pn{a)v = 0 is of dimension p. If we write Xn,k, 
1 < k < m, for the different zeros of Pn, and Ik is the multiplicity of Xn,k, then 
the matrices 

Tfc = , , \ ^ {Xn,k)Qn{Xn k): 1 < fc < 771 

(det(P„(t)))i‘'“Tx„,fc) 

are positive semidefinite of rank Ik, and with them the matrix quadrature 
formula 

/ m 

F(t)dVk(t) = ^ p(x„,fe)r„,fc 
k=l 

holds for every matrix polynomial P of degree at most 2n — 1. 

If we assume that An —> A, Bn B where A is non-singular, then 


Pn{z)P-l,{z)An^ 


dWA,B{t) 

z — t 


locally uniformly outside the cluster set of the zeros, where WA,B{t) is the mea¬ 
sure matrix of orthogonality for the sequence of matrix orthogonal polynomials 
Sn with recurrence coefficients A, B for all n, i.e., which satisfy the three-term 
recurrence 

tSn{t) = A*Sn+l{t) + BSn{t) + ASn-l{t). 

The distribution of the zeros themselves will be 1 /A^-times the trace of the ma¬ 
trix measure of orthogonality for another sequence of matrix orthogonal poly¬ 
nomials Rn satisfying 


tRn{t) — ARn+l{t) + BRn{t) + ARn-l{t), 77 > 2, 


with appropriate modifications for n < 2. 


Families of orthogonal polynomials 

If the matrix of orthogonality is diagonal (or similar to a diagonal matrix) with 
diagonal entries pi, then the orthogonal matrix polynomials are also diagonal 
with 7-th entry equal to Pn(pi), the 77 -th orthogonal polynomial with respect 
to Pi- Many matrix orthogonal polynomials in the literature can be reduced to 
this scalar case. Recently however, some remarkably rich non-reducible families 
have been obtained by A. Duran and F. Griinbaum (see [25] and the references 
therein), which may play the role of the classical orthogonal polynomials in 
higher dimension. They found families of matrix orthogonal polynomials that 
satisfy second order (matrix) differential equations just like the classical orthog¬ 
onal polynomials. Their starting point was a symmetry property between the 
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orthogonality measure matrix and a second order differential operator. They 
worked out several explicit examples. Here is one of them: N = 2, the measure 
matrix (more precisely its density) is 


H{t) := e-*' 


1 + at^ \ 

at^ 1 y ’ 


t G H, 


where a G C \ {0} is a free parameter. The corresponding P„(t) satisfies 




—2t Aat \ 
0 -2t ) 


Pn{t) 


-4 2a \ 

0 0 ) 




There is an explicit Rodrigues’ type representation for the polynomials them¬ 
selves, and the three-term recurrence (16.1) holds with Bn = 0, 


An+l 


n + 1 / 

2 V 0 


<l'^n+2ln+l 

Inhn+l 


where 



Connection with higher order scalar recurrence 

Matrix orthogonality is closely connected to {2N -|- l)-terni recurrences for 
scalar polynomials. To describe this we need the following operators on poly¬ 
nomials p: if p{t) = ^ then 

PN.miP^ — ^ ^ O^sN+mil t 

s 

i.e., from a polynomial the operator RN,m takes those powers where the expo¬ 
nent is congruent to m modulo N, removes the common factor and changes 
to t. 

Now suppose that {pn}?t=o is a sequence of scalar polynomials of correspond¬ 
ing degree n = 0,1,..., and suppose that this sequence satisfies a (2A^-|- l)-term 
recurrence relation 

N 

t ~ ^npPnit) T ^ ^ (Cyi.fcPn —fc (Q T Cn-\-k,kPn+k{l^\ 

fe=l 




V. Totik 


118 


where c„^o is real, c„^jv 7^ 0 (and Pkit) = 0 for fc < 0). Then 


Puit) = 


( RN,o{PnN) 
RN,oiPnN+l) 

V RN,o{PnN+N-l) 


RN,N-liPnN) 
RN,N-liPnN + 1 ) 

RN,N-l{PnN + N — 1) ) 


is a sequence of matrix orthogonal polynomials with respect to a positive definite 
measure matrix. Conversely, if Pn = {Pn,m,j)m~j=o i® ^ sequence of orthonormal 
matrix polynomials, then the scalar polynomials 


N-l 

PnN+m{t) = ^ PPn,m,j{t^), 0<m< N, n = 0,1,2,..., 
j=0 


satisfy a {2N + l)-recurrence relation of the above form. 
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